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jL . Abstract 

i-C ' Abstract 

We examine the generation of baryon asymmetry at the weak scale from a primordial 
K> ' lepton asymmetry. If the electroweak phase transition is first order, partial reflection 

C^ ' of tau leptons off the bubble walls and the resulting hypercharge transport into the 

symmetric phase, gives rise to a baryon number consistent with observation. 



Introduction 

Electroweak theory contains the means for anomalous violation of fermion number via non- 
perturbative effects, in the form of instantons at zero temperature |jl| and sphalerons at 
high temperature [^ . The vacuum structure of electroweak theory admits both quantum 
and thermal tunneling between inequivalent vacua by such configurations. Because of 
this, any theory in which i? — L = is in danger of having its baryon number erased by 
sphalerons. 

One way around this problem is to suppose that the baryon asymmetry was made at 
the electroweak scale [Q Q . If this is the case the sphaleron transitions must have fallen 
out of equilibrium whilst there was still some chemical potential for non-zero baryon 
number. If the CP violation is provided by the higgs sector, this implies that the phase 
transition was at least mildly first order [^, Q. 

The alternative possibility which will concern us here, is that the asymmetry was 
encoded in a conserved global quantum number such as B — L number 0. This was 
generated (or set by hand) primordially, and after possible processing at the electroweak 
phase transition, was converted into the presently observed baryon number. Recent work 
@, H has suggested that the three 

L^ - B/3 

numbers may play this role {i is a generation index). These could have been produced by 
generational (and CP violating) differences in the out of equilibrium decays of lepto-quarks 



P, |I3]. The baryogenesis could have occurred after the electro- weak phase transition as 
a results of mass effects. For this to be possible it is necessary that sphalerons remained 
in equilibrium for some time implying a second-order transition. Campbell, Davidson, 



Ellis, and Olive |11] have taken this idea further by demonstrating how to avoid the 
im? /T"^ suppression through lepton- violating interactions. Provided that some, but not all, 
lepton flavours are violated by AL 7^ interactions in equilibrium, B may be regenerated 
without the m^/T^ lepton mass effects. Instead B depends only on the initial asymmetry 
of the non-equilibrating modes. For example, if lepton-violating interactions equilibrate 
Li — 5/3 and L2 — -B/3 to zero, then the non-zero L3 — B /?> biases baryon production 
via sphalerons: 

B ocB - L = B/3 - L3 ^ . 

The above example is interesting for two reasons. Firstly it can operate even ii B — L = 
0. Secondly there is the novelty that two of the Sakharov conditions (C/CP violation and 
thermal non-equilibrium) are important at early times, whilst the third (baryon number 
violation) is important at 'late' times (where 'late' means at the weak scale; after ^ 10~^° 
seconds). In the light of this, we shall discuss an alternative way in which these quantum 
numbers can furnish a net baryon number. 

Again we rely on some unspecified mechanism to generate a non-zero L^—B/S number. 
The values of the other two numbers are unimportant, and it could be that B — L = 0. 
The phase transition is assumed to be robustly first order, as argued for in the case of the 



Standard Model in refs.H, 121, and as demonstrated in extensions of it in ref. 13 . Once 



the temperature drops below the critical temperature, bubbles of true vacuum begin to 



nucleate and expand into the symmetric phase. As the walls advance through the plasma, 
a fraction (proportional to their physical mass) of the r-leptons are reflected back into 
the symmetric phase. Since the L3 — 5/3 number is non-zero, more taus than anti-taus 
are reflected, and a small hypercharge builds up on the outside of the bubbles with an 
opposing hypercharge on the inside. On the outside sphaleron transitions convert some of 
the lepton number into baryons, whilst on the inside sphaleron transitions are suppressed 
by a large exponential factor, resulting in a net number of baryons being produced. 

Since the reflection is only suppressed by the factor Mr/T one might expect this 
mechanism to be much more efficient than the equilibrium process described in refs.[§, ^. 
We shall see that in fact the two mechanisms give remarkably similar baryon asymmetries. 
This is due to an additional suppression due to the low rate of conversion of lepton into 
baryon number compared to the typical wall velocity. The resulting baryon number is 
proportional to the initial r-lepton asymmetry p^-. We stress that we do not make any 
assumption about enhanced (or maximal) CP violation in the dynamics of the phase 
transition (which plays no role in our mechanism 0). 

At high temperature or density, the fundamental excitations do not coincide with the 
elementary particles. As a result, perturbation theory in terms of the bare fields does not 
properly describe effects essentially due to the ambient plasma. In order to proceed in 
a valid perturbative calculation, we first determine the quasiparticle modes that exist in 
a relativistic plasma at high temperature. Using these modes, we compute the refiection 
coefficients and then the lepton fiux refiected off the expanding bubble wall. Because the 
reflected flux oc Mi/T, where M; is the lepton flavour mass and T ~ 100 GeV is the critical 
temperature of the phase transition, the r-lepton flux will predominate. The thermal 
scattering length of the leptonic quasiparticles will then be computed, to give an indication 
of how long the reflected leptons remain in the unbroken phase before being absorbed by 
the advancing front of broken phase. During its time in the unbroken phase, the r-lepton 
flux biases anomalous baryon violation, which, it is assumed, are in equilibrium outside 
the bubble (the region of = 0). The sphaleron processes thus produce baryons, which 
are swallowed into the broken phase and survive to the present. We then estimate the 
generated baryon number and flnd that pb ~ 10~^° can be accommodated in this scenario, 
provided that the primordial r-asymmetry is Pr/s ~ 5 ■ 10~^, for a wall velocity u ^ 0.1. 
Scaling by the lepton number density, 

^^i^:^ ~ 0.005 . 

ni + ni 
Our primary sources are the papers by Nelson, Kaplan, and Cohen [Q and by Farrar and 



Shaposhnikov [|12[, denoted NKC and FS respectively. The starting point of our analysis 



is the rate equation describing the approach of baryon number to equilibrium |[T^ ; 

Tb of 

^CP violation is of course evident in the initial condition of a primordial L-asymmetry; and CP 
violation in the particle dynamics would have been required for its generation at the GUT scale, as 
required by the Sakharov conditions. 



The partial derivative of the free energy is taken with all conserved quantum numbers 
held fixed, and as we will compute later, is simply proportional to the hypercharge; 

— = e^ (2) 

where the parameter ^ depends on the model under consideration and will be determined 
explicitly later. Integrating the rate equation ahead of the advancing bubble wall, 

r^ f OF ^Tb 1" , ^. 

Pb = —^^ cit - — = — T—r I at Py[z — ut) 

— CO 



where fy is the reflected hypercharge flux and tt the thermal transport time (~ thermal 
scattering length). Scaling by the entropy density and recalling the rate for anomalous 
baryon violation in the symmetric phase (for three generations) [jT5|, |1^, 



"b = 3Ka%T'^ , (4) 



we arrive at our final expression for the observed baryon number, 

Pb 3a w^^ f fYTTT\ 



s u 



(5) 



We consider the ranges 0.1 < k < 1 [T^] and 0.1 < m < 1 [jT^, and are left to compute /y, 
^, and tti which we do in the following sections. For the remainder of this discussion, we 
shall assume for concreteness that i? — L = 0, although this is not necessarily the case as 
stated previously. 

Reflected Hypercharge flux (/y) 

Before we can calculate fy we need to determine the leptonic modes in the thermal plasma 
(which is characterised by the four-velocity -u") and obtain the quasiparticle solutions for 
left and right chiralities, which have different interactions and hence develop different 
thermal masses. The first part of this section section summarises basic results; more 



complete treatments may be found in Weldon and Lebedev [18 
The lepton thermal self-energy at one loop is given by 

S(P) = ^AJ-0^^D,^{k)rS{k + P)Y + ^BJ-0^^D{k)S{k + P) (6) 

= -a{uj,p)^-b{uj,p)]i, (7) 
where a and b are functions of the Lorentz invariants 

iu = P-u (8) 



p ^ J{P.uf-P\ (9) 



such that uj"^ — p^ = P^. The functions a and h are computed in |jT8|: 






io I UJ +p 

2p \^-p, 



b{uj,p) 



p 



u 1 

— + 



u 



p 2 \p^ 



- 1 l0£ 



UJ +p 



1 T2 



(10) 

(11) 

(12) 



The constants A and B depend on the lepton chirahty and the model under consideration. 
For the minimal Standard Model and the two-doublet model, one can easily see that 



A, 



h' + y 



AR = g 



a 



Bl 
Br 



9r 



(13) 



where the Yukawa contributions are identical for both models, since in either case, the 
lepton couples to only one scalar doublet. This yields the following thermal lepton masses: 

1 



nl 



—- (2M^ + M| + M2 T^ ^ (0.209 T)2 



nl = -L(2tan2^wM^ + M2)T2^ (0.127 r)^ 



(14) 
(15) 



Note that Ql > ^r] this is true for all leptons /, provided that Mi < 116 GeV. Also, this 
result is gauge invariant. The self-energy (^ leads to the following fermion propagator: 



S{P) = [{l + a)^+b: 
where the denominator is 



[{l + a)f^+bji]/Z 



(16) 



Z{uj,p) = (1 + ayP' + 2{l + a)bP-u + b' = [uj{l + a) + b]' - [p{l + a)Y ■ (17) 

The poles of the propagator occur at the zeros of Z{ijj,p), 

uj{l + a) + b = ±p{l + a) , (18) 

which includes both positive- and negative-energy solutions (note that (1 + a) is even, 
and b odd, under oj -^ —uj); given a positive-energy solution u^p), the corresponding 
negative-energy solution is —uj{p). So it suffices to find only the positive-energy solutions. 
The dispersion relation is given by 



UTP 



-a{uj,p){uJTp) -b{uj,p) 
fi2 



P 



-^K-?'-^, 



(19) 



Each chirality has two distinct modes, shown in Figure 1, which we label normal {uj^{p)) 
and abnormal {uj_{p)). The abnormal mode is actually unstable for p ^ Q [|1^. We 



plot the dispersion relations in Figure 2 for left and right chiralities of the r-lepton. The 
dispersion relations may be approximated for small and large p as 



^+(p) 



1 






i p. 
3 \n 






1 p 
3 \n 



1 + 2 exp -2 



: p ^ Q, 




: p ^ ^ 






: p ^ Q, 




: p ^ il 



(20) 



f211 



In the broken phase, the Dirac operator is 



y^ 

^L 



y^ 



M 



Mt y% + y\ 



:22i 



where the thermal piece includes a contribution due to mass corrections in the broken 
theory: 

(23) 



'L,R 



^L,R + ^^L,R 



We then obtain the dispersion relation from 



det ( ^° 



Mt 



^b 

^L 



M 



0. 



(24) 



The Lagrangian may be made linear for p <^ ^l,r, 



c 



eflf 



1 



2iL^ (do--a-d + iQL]L + 2iR^ ido + -ff-d + iQR]R 



1 



+ L^MR + R^M^L 
yielding the dispersion relation 

^L + ^R 



(25) 



uj{p) 



± 



M2 



'Vti 



-^R^\P\' 

2 3 , 



(26) 



where the first ± is for L, R and the second for normal or abnormal, respectively; the 
solutions uj{p) are shown in Figures 3 and 4. Notice that the left abnormal and right nor- 
mal lines do not intersect (as they do in the symmetric phase), but instead are separated 
by an energy interval Auj = Mj. about the point cjq = (^l + ^r)/2. This interval is the 
region of total reflection, as we will see below when we consider the scattering of fermions 
off the bubble wall. 



In order to calculate the reflection coefficients, we start from the full Dirac equation. 



yl + El + sy^ 



^% 



M 

^% + syi 



L 
R 



0, 



(27) 



where Ai (x) is the position-dependent mass and embodies the details of the bubble wall 
profile; in the broken phase S*^ = S*^ + 5S and Ai = M, while in the symmetric phase 
STi = = Ai. It is because the left and right chiralities interact differently with the bubble 
wall that there exists the possibility of seperating C- and CP-odd reflecting currents (weak 
hyper charge for instance). 

As a first approximation, we set S^ ^ S*^ or 5S ^ 0, which is plausible since for 
leptons M/T <^ 1; (5S includes mass corrections in the one- loop graphs via the propagator 
\^ + M]~^, where the momentum integrals are dominated by the region p ^ T. Again we 
may linearise the Dirac equation at low momenta, p <^ Q: 



ijj 



(l + aL + Z^i) +ia-(9(l + aL) M i(-i_q 




This may be written as 



\ 9^ 



iUR^i , (29) 



where 

R = f^'-'r^" . ° .-A (31) 




* = ^'M D I (32) 

OIL = -^'-^{l-3u-2u')il-u) ; PL = ---^il+ufil-u) (33) 

3 (jJ 

aR = -^'-^il + 3u-2u'){l + u) ; Pj, = -1^(1 - u)\l + u) .(34) 

3 cu"^ ^ ^ 3 cj"^ 

A Lorentz transformation from the fluid frame to the wall frame (with velocity u) has 
been performed in the small u limit; the -u-dependence of a and (3 reflect the spinor 
transformation. Because we are working in the wall frame, the energy and momentum 
parallel to the wall are conserved: i-^ = u, i^- = p\\. The reflection coefficient depends 

strongly only on p^ (taken to be Pz), and so we set p\\ = 0: 




iUR"^ , (35) 



which decomposes into 



d 


' i'l 


dz 


,^3 


d 


'^. 


dz 


, V^2 



Vl + tVL/ V l+CVfl / 



1 + Qfl ) \1 + Ql, 




(36) 

(37) 



This describes reflection and transmission perpendicular to the wall. In analogy to the 
method of NKC, we write the solution as a path-ordered integral, 

where 



n{z) = Vexpi I dx{ ^ /+"^, { /iVt"«^ X I (40) 



20 



MUr^ 



We denote the path-ordered exponentials by 

2 

n{z) = Vexpi fdxQ{x) (42) 

-20 

Z 

U{z) = Vexpi fdxQ{x). (43) 



-20 

Note that fi and fi satisfy the differential equations ^ = iQ and ^ = iQ- ^ describes 

L ^ R reflection and R ^' R transmission, while f2 describes R -^ L reflection and 

L ^ L transmission. As a first approximation, we take the wall profile to be described 

by a step-function: 

i M : z>0 ,,,, 

-^ = 1 n n • (44) 

Then in the symmetric phase, 

^ ^ \ exp -2CU (i^^fg^U / V exp-ip%zj 

while in the broken phase, Q may be diagonalised as 



(45) 



n{z>0) = D-'i ^"P^^^" °., D, (46) 

'0 exp —tpj^z I 



where the momenta in the broken phase are 



Pl - Pr 



{ pi+p% y M^ 



^^^•^" 2 ^\[ 2 ) {l + a,){l + an 



(47) 



Then 



D 



^1 



(z>0) 



exp ipi^z 



,^3 J' ' \ 

with the diagonahsation matrix given by 

^U+VB) 

X \l+aL) 

where 

_ pI+Pr 

^ — z 



D 





exp —ip\z 






D 



^3 



(0) 



B 



u^ 



M2 



{I + aL){l + diR) 



(48) 

(49) 

(50) 
(51) 



Here X and Y are just normahsation constants. The reflection coefficient conies from the 
condition 

^ ^^ d{ \\ , (52) 



yielding 



n 



L^R 



\R\ 







L~^R 



l+aR 



1 + aL 



UJ 



uJ+VB 



(53) 



Similarly, the right-to-left reflection coefficient may be found from the Dirac equation (|37|) 
for iIj2 and ip^, with the result 



7^ 



R^L 



l + ttL 

l + aR 


oj-VB 


uJ+y/B 



Note that TZl^r ~ TIr^l ~ 1 if VB is imaginary; that is, if 

M2 



0J'< 



(54) 



(55) 



;i + aL)(l + aR) ' 

we get total reflection. In physical terms, an imaginary contribution to the momentum 
(^) implies an evanescent (decaying exponential) transmission amplitude in the broken 
phase. 

Although the step-function wall proffie may produce spurious effects due to its discon- 
tinuity, Farrar and Shaposhnikov considered the smooth proffie 

M2 



M' 



1 + e- 



(56) 



and found the reflection coefficients 

TZl^r = 

TZr^l = 




sinh - (uj 



sinh^(cJ+ VB 



sinh - (uj 



sinh^iuJ+VB 



(57) 



Notice that total reflection again appears for ufi < M^/(l + (5l)(1 + au). We will take 
this as a general condition for total reflection (like the condition u < M for the zero- 
temperature case). The region of total reflection is shown in Figure 5. As expected by 
comparison to the dispersion relations in the broken phase (Figure 4), total reflection 
occurs in the energy interval of width Aa; ~ M about uq ~ {Ql + ^_r)/2; the deviation 
results from the dependence of the reflection coefficients on the wall velocity u. 

Using this result we are now able to compute the hypercharge flux reflected off the bubble 
wall. We neglect quarks, since we assume that the primordial baryon asymmetry vanishes 
and that CP violation in the CKM matrix is inadequate to produce the observed baryon 
number. We need to calculate 



+ (/i-contribution) + (e-contribution) . (58) 



Because M^ ^ M^, Mg, the interactions of the muon and electron families with the bubble 
wall are negligible compared to those of the tau family, and hence their contributions to 
the reflected flux may be ignored. In the following, /^-^ denotes the tl particle flux in 
the fluid frame, while f^ denotes the tl number density distributions in the wall frame, 
for the symmetric and broken phases; and similarly for the other species. To flnd fj, 
the particle flux in the thermal (fluid) frame, we compute the flux 7/^ in the wall frame 
(7 = 1/Vl — u'^, where u is the wall velocity): 

j-^ [/^(A;^,M-^i?^L(^)] (59) 

^^ [mkL,kT)-nL-.RH] . (60) 

The integrals are taken over particle momenta in the wall frame and by CPT and Lorentz 
invariance, TZi^ji = TZl->r- There is no term for reflection off the symmetric phase 
because the dispersion curves for left-abnormal and right-normal modes intersect in the 
symmetric phase; see Figure 2. Also there is no term for particle transmission from the 
broken phase since it is assumed that both the broken and symmetric phases have zero 
net hypercharge, so that the transmitted hypercharge flux is zero when summed over 
all particle species. Thus reflection yields the only substantial hypercharge flux, and 
in particular, the contribution from total reflection dominates. The difference of these 
integrals may be written as 

1 r rl^k 
fr, -/-« = - y (^ IfR C^n-L - n,^n) + ifn - /!) -^l^r] ■ (61) 

The flrst term is much smaller (by a factor of ~ 10) than the second since the difference 
in reflection coefficients results from the difference in the left and right thermal masses. 
We shall therefore omit it in the remainder of this discussion for the sake of clarity. 
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The above integral is dominated by the region of total reflection, 7^ ^ 1, which (as we 
have seen above) occurs for i? < or 



UJ'< 



and the hypercharge flux is therefore governed by the r-lepton; 



(62) 



/y ~ 



' r) kJtl Jtr) {Jtj{ It 



7 J (27r)3 
1 f d^k 

7 



B<0 






The flux distributions are taken in the wall frame; 



f^ 



dujj 



Up 



b 



u, - uP; 



Tf^j 



(63) 



(64) 



'' dk, 
where the group momentum and fermion particle distribution (in the thermal frame) are 






Substituting into eq.(|63D gives 
1 f d^k f 1 dijR 



h 



7 J (27r)= 



[np {-i[ujR - mP|] - ^Itr)- riF {i[ujr - uP^] + fim) 



(65) 
(66) 



2 dk, 

+ -Qj^ [np hiuJL - uP[] - firj- np (7K - uP[] + /i^^)] 



(67) 



As a conservative estimate, we truncate the region of phase space where -B < to the 
region where k\\ ^ |(^l — ^b) and dujLRJdkw ~ (see FS for a discussion on this point). 
The k\\- and /c^-integrals can now be separated, giving the result 

fy ~ -^lTZ7,^^i^L-^Rf\-:^^nF{l[uJQ-uP|^]T^^rn) 



167r7 



+ Anp {i[uJo - uPI]t ^^rS\ 



(6^ 



where the region of total reflection is centered about uq with spread Acu. To lowest order 
in the chemical potentials /i^-^^T^^/T, 



/^np (7 cjo — uP^ 



T/ij 



) = np[ 






f^jj 



) - ^F (7 ^0 - uPj + l^j 



cosh 



■y{u}0-uPf) 



(69) 



11 



and in the limit of small u and k\\ , 

P^^^(l + 2«/3), (70) 

hence 

iA„. (. I.. - UPS ^ ,,) . A„. (-, [.. - .Pn ^ ..) 4 ■ J:2l^^l^ ■ 

(71) 
Noticing that Pr = {Pm + '^f^TL)T'^/Q (refer to the following section), we find that fy 
is simply proportional to the primordial r-lepton asymmetry. Our final result for the 
reflected hypercharge flux is then 

167r T3 i + cosh(f(l-M/3)) ■ 



We believe this expression to be a conservative lower bound, since we have truncated the 
flux integral (|68|) in the region of total reflection and we have neglected the contribution 

of fc|| ~ Uq. 

Partial Derivative of the Free Energy (dF/dB) 

We now compute dF/dB, the partial derivative of the free energy with all conserved 
quantum numbers held flxed. We cannot simply adopt the result of NKC, since their 
analysis implies not only B = L = 0, but also Bj = Lj = for the individual generations. 
Since we are interested in the case of generational differences in lepton asymmetries, we 
redo the analysis. 

dF/dB depends on the interactions and species in equilibrium. We assume that on 
the timescale tt, SU{3)c ® SU{2)l ® U{\)y interactions are in equilibrium (including 
quark mixing and light-fermion Yukawa interactions, which are assumed to be out of 
equilibrium in the NKC analysis), and that only the anomalous B + L- violation is not 
in equilibrium. The (approximately) conserved quantum numbers and their associated 
chemical potentials are then 



-L, 


^ /ij 


B 


^ Pb 


Y/2 


^ I^Y 


Ts 


^ Pt 



This implies that on the timescale of interest, Q, Lj, and B ± L are also conserved. 
In contrast NKC had Bi = B2 and B3 separately conserved, and Lj = L/3; when the 
constraints B = = L and Bi = B2 = were imposed, one obtained Bj = = Lj. In our 
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case we have Lj conserved and Bj = B/3, and the constraints we impose are B = = L, 
so that Bj = although Lj 7^ is allowed. The net particle number densities are 

rl , ,1 1 1 1 T^ 

- (^1 + /i2 + /is) + g/iB + qPy + ^Pt 



PtL 

PbL 

PtR 

PbR 
PCL 
PCR 

PeL 
Pt^L 
pTL 
pen 



1 , ,1 1 1 ■ 

- [Pi + P2 + P3) + -^Pb + -Py - -Pt 



6 



1 12 

- (/il + P2+ Ps) + -PB + -pY 

1 , ,11 

- [Pl + P2+ P3> + -pB - -pY 

puL = PtL ; 

Pub ~ PtR 1 



Tl 

6 

Tl 

6 

Psi = PdL = PbL 

PSR = PdR = PbR 



-Pl-^PY--PTJ g 






1 1 

2^^ ~ 2' 
1 1 

2^^ " 2' 
1 1 

^2 



-/i2 - 2Py - 2/^T^ 6 



Pu, = 


( 


1 1 \ 
-^PY + -/iT j 


T2 
6 


Pu^ = 


( 

-P2- 


1 1 \ 

- 2/^y + -^Pt] 


T2 
6 


Pur = 


( 
{-P3- 


1 1 \ 

2/^1- + IjPT ^ 


T2 
6 



-yUl - yUyJ 



P0+ = n (/iy + /Xt) 



6 

rn2 



P/^fl = (-P2 - py) 



6 



'T-R 



;-p3 - py) 



T2 



n {iiy - [It) 



rn2 



Pw+ 



^P7 






In the above, n denotes the number of scalar doublets in equilibrium. Then 

PY I , ^ 2 1 



- (peR + Pm + Prfl) + - (P0+ + P<^o) 

Q 

- (pi + P2 + Pa) + 2pB + (10 + n)fiY 



11 
6 



PT3 = (10 + n)fiT 



6 



Pb 



PL 



3- 2 (PtL +P6l +Ptfl+Pbfl) 

'4 

- (Pl + P2 + Pa) + 4/iB + 2/iy 



6 



PcL + Pl^L + Ptl + Pue + Pl^M + P^r + Pcfl + P^fl + pTR 

T2 



[-3 (/ii + H2 + Ps) - 6py] 



6 
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rp2 rp2 

PL, = (-3/ij - 2/iy) — = [-3/ij + (/ii + /i2 + /is)] Y • ( ''3) 

Imposing the conditions B = = L, we find 



giving 







flY 


= 


6pB 


/il 


+ ^2 + ^3 


= 


-12pB , 


9F 


= l^B - 


/ r' 


PY 


_ ^PY 

mO mO 



(74) 



Thermal Transport Time (tt) 



(iViv) = (iy/^) , (75) 



After rebounding off tlie bubble wall, the reflected lepton flux travels into the symmetric 
phase until the advancing front of broken phase captures it, during which time the cor- 
responding hypercharge current biases baryon production via anomalous processes. We 
now estimate the thermal transport time tt, defined as the average time that a reflected 
r-lepton spends in the plasma prior to absorption by the bubble of broken phase. Consider 
diffusion away from the wall of a particle with velocity v and mean free path /. Capture 
occurs when the wall intercepts the randomly walking particle (A^ here is the number of 
collisions) : 

fvTr\ 

UTt '■-■■'• 

or 

We therefore need to calculate the thermal average 

(M = (t^\ , (77) 

where the subscripts 1 and 2 refer to thermal averages taken over the rebounding particles 
and particles in the plasma, respectively. Here n is the particle number density and a the 
thermally averaged cross section. The thermal transport time may then be written as 

-1 _^ U^9l92 r d^Pld^P2 El CTl2iPl,P2) ,„^x 

where the sum is taken over all interactions of particle 1 with the heat bath and Qj counts 
the spin degrees of freedom. Since the leptons only interact electroweakly, it is reasonable 
to approximate the above expression by considering only the leading contributions from 
tree-level scattering on mass shell. Numerical calculation confirms that other contribu- 
tions are indeed less significant. In this case, the transition probabilities sum to 

Em' - -iri^'w-^l) (79) 

sm t^w ^ ^ 
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+ 4e^ {Q^^ + Q 



2' 



(81) 
(82) 



where we respectively list the contributions of W, Z, and 7 scattering in the plasma (higgs 
scattering is negligible since it is suppressed by the small Yukawa couplings). Qj is the 
thermal mass of particle species j, and we have set Q^ ~ Qw ~ 0.5T and fij,^ ^ Q^ ^ 0.2T. 
The thermal transport time is then given by 



X 



100 



Tt = 



u^T u^T 



(83) 



Note that, computing (f ) = {nv)/{n), we find that v ^ 1 for most leptons in the plasma, 
which confirms that the reflected lepton current is quickly thermalised. 

Discussion 



Taking the expression (|^) for the final baryon number and substituting in eqs.(|72]), (|7^, 
and (|8^) for the reflected hypercharge flux, the partial derivative of the free energy with 
respect to baryon number, and the thermal transport time, we obtain 
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(84) 



where the various parameters have been scaled by their typical values. We believe eq.( 
to be a conservative lower bound on the effect of this mechanism since we have underes- 
timated the flux integrals. By requiring the observed baryon number of ps ~ 10~^° to be 
generated in this manner, we obtain an estimate for the primordial lepton asymmetry. 



Pt 



5-10" 



for the typical values of parameters (k = 1, u = 0.1, x = 100). This corresponds to 

ni - nj 



ni + nj 



0.005 . 



^5) 



(86) 



It is of interest to compare this constraint with that obtained from equilibrium scenar- 
ios for lepton-to-baryon conversion. As a generic example, we consider the analysis of 
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Kuzmin, Rubakov, and Shaposhnikov PJ. Taking the requisite large Higgs mass to be 
Mh ~ 100 GeV, and consequently the sphaleron freeze-out temperature (approximately 
the critical temperature) to be T^, ~ 150 GeV, we estimate the generated baryon asym- 
metry as 

P^ - ^ "^' ^^) ^ -4.4 . 10-^ ■ (P-) . (87) 



s 137r2 T2 

Hence the observed baryon asymmetry may be accounted for in this scheme for a primor- 
dial lepton asymmetry of pr/ s ~ 2 ■ 10~^, which is remarkably similar to the value found 
above despite the ml/T"^ suppression. This is not surprising since, as we have emphasised, 
the assumption that baryon violation is in equilibrium (on the timescale of interest) is an 
overly optimistic one. 

A remark is in order about the behaviour of eq.(0) as m ^ 0. In contrast to NKC, 
our result diverges unashamedly at low wall velocity as ~ 1/u^. To see why this is 
acceptable, consider the extreme case of a stationary wall. In this case the r-leptons 
simply diffuse away from the wall, finally establishing equilibrium once the hypercharge 
profile is constant (and different) on either side. In this sense the wall is behaving like 
a semi-permeable membrane since it forces different partial pressures of leptons on the 
outside and inside of the bubble. Intuitively, one expects the partial pressure of leptons in 
the symmetric phase to create baryons via sphaleron transitions until it is fully depleted. 
However in this calculation we have made an approximation by setting the chemical 
potentials to be constant. This means an unlimited supply of lepton number - hence the 
divergence. Our approximation breaks down once the wall velocity is slow enough so that 
the lepton numbers are depleted significantly. This condition can be gauged by the ratio 
of generated density of baryons to initial density of leptons. This is a small number for 
u ^ 0.001 which gives a sufficient range of values for first order phase transitions (where 

u ^ 0.01 is more reasonable) Jl^ . Clearly, in this scenario, the slower the wall velocity 
the more efficient the mechanism. 

To summarise, we have considered the generation of baryons from a primordial r-lepton 
asymmetry. There are two ways in which baryogenesis may occur which give remarkably 
similar results. The first is already well known in the literature, and assumes that the 
thermal plasma maintains equilibrium during the electro-weak phase transition. Anoma- 
lous processes may convert generational lepton asymmetries into baryon number, whose 
final value depends on the primordial asymmetry. The effectiveness of this scenario is de- 
termined by the suppression factor Mf{T^)/T'^ at sphaleron freeze-out. We have analysed 
a second mechanism, which would operate during a first-order phase transition. Refiection 
of r-leptons off the phase separation boundary may radiate a net hypercharge flux, which 
then triggers baryon production as described by the rate equation (|l]). The effectiveness 
of this charge transport mechanism is determined by the strength of the lepton Yukawa 
interactions with the bubble wall, through the factor Mi/T. This mechanism works most 
efficiently for slow wall velocities. Although a comparison between the two depends on 
the choice of parameter values, there is a clear trade-off between opposing tendencies: 
anomalous baryon violation in equilibrium generates greater net 5, but risks suppression 
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by lepton mass effects. We have made several assumptions in deriving the final baryon 
number, which we summarise: 

• the reflected hypercharge flux is dominated by the contribution due to total reflec- 
tion; 

• the linearised (low-momentum) Lagrangian is valid in the region of total reflection; 

• the wall velocity is non-relativistic (implying that the hypercharge current rapidly 
thermalises); 

• the flux integrals may be approximated by k\\ ^ |(^l — ^r)- 

The last assumption together with the parameters k, u and x are the greatest unknowns. 
For typical values, we have found that the observed baryon number of pb ~ 10~^° may 
be generated by our mechanism, if the primordial lepton asymmetry is Pr/ s ~ (9(10~^). 
Finally, we note that scalar leptons in supersymmetry, may play a similar role to the 
r-lepton above. After finding the quasiparticle modes and dispersion relations by diag- 
onalising the mass matrices in both the symmetric and broken phases, one may analyse 
the scattering of sleptons off the bubble wall in the fashion described above. We expect 
that a region of total reflection of width Auo ~ M^ would again yield an M^/T-suppressed 
contribution to the reflected hypercharge flux, in the manner demonstrated. This is in 
contrast to the equilibrium case examined in ref. 0, where the minimal supersymmetric 
extension of the Standard Model produces a much greater enhancement due to the large 
mass splittings of right and left sleptons. 
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FIGURE CAPTIONS 

Figure 1 : The dispersion relations for normal and abnormal plasma modes in the symmetric 
phase. 

Figure 2 : The dispersion relations for left and right chiralities of the r-lepton in the symmetric 
phase of the Standard Model and two-doublet model. 

Figure 3 : The dispersion relations, at low momenta in the broken phase, for normal and 
abnormal modes of left and right chiralities of the r-lepton in the Standard Model 
and two-doublet model. 

Figure 4 : Close-up of the previous figure, magnifying the separation between the left abnormal 
and right normal lines. 

Figure 5 : The region of total reflection (bounded by the curves plotted) for the r-lepton in 
the Standard Model and two-doublet model, as a function of the wall velocity; the 
temperature is taken to be T = 100 GeV. 
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